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Abstract. We study a moving boundary problem describing the growth 
of nonnecrotic tumors in different regimes of vascularisation. This model 
consists of two decoupled Dirichlet problem, one for the rate at which 
nutrient is added to the tumor domain and one for the pressure inside 
the tumor. These variables are coupled by a relation which describes the 
dynamic of the boundary. By re-expressing the problem as an abstract 
evolution equation, we prove local well-posedness in the small Holder 
spaces context. Further on, we use the principle of linearised stability 
to characterise the stability properties of the unique radially symmetric 
equilibrium of the problem. 



1. Introduction 

The study of tumor growth models is a very current topic in mathematics. 
During the last four decades an increasing number of mathematical models 
have been proposed to describe the growth of solid tumors (see [3,7,9,14] 
and the literature therein). There is a three level approach in modeling the 
complex phenomena influencing and describing the processes inside a tumor. 
Models at sub-cellular level take into consideration that the evolution of a 
cell is determined by the genes in its nucleus, at cellular level they model 
cell-cell interaction and at macroscopic level, when the tumor is considered to 
consist of three zones: an external proliferating zone near high concentration 
of nutrient, an intermediate layer and an internal zone consisting of necrotic 
cells only. Very often models combine aspects from different scales. There are 
also, a large variety of different types of models: biological models, consisting 
of coupled ODE systems where the variables correspond to some biological 
properties of an entire population; mechanical models yield to determine the 
cell movement based on physical forces; the discrete models handle single- 
cell scale phenomena and the effects are then examined at macroscopic scale 
and moving boundary models when the macroscopic description of biological 
tissues is obtained from continuum mechanics or microscopic description at 
cellular level. 
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In this paper we deal with a moving boundary problem, which is obtained 
by combining aspects from the cellular and macroscopic scale, and possesses 
also characteristics of the mechanical model (Darcy's law). Cristini et al. 
obtained in [6], using algebraic manipulations, a new mathematical formu- 
lation of an existing model (see [5, 14, 18]), which describes the evolution of 
nonnecrotic tumors in all regimes of vascularisation. This new formulation 
has the advantage of considering different intrinsic-time and length-scales 
related to the evolution of the tumor and, by incorporating them in the 
modeling, provides a model describing both vascular and avascular tumor: 

in n{t), t > 0, 
in n{t), t > 0, 
on dn{t), t > 0, 

I |2 

AG^-^ ondn{t), t>0, 
on d^l{t), t > 0, 

(1.1) 

Hereby J^o is the initial state of the tumor, V is the normal velocity of the 
tumor boundary, HQQ(^t) the curvature of d^{t), and the constants A and G 
have biological meaning, namely G is the rate of mitosis (cell proliferation) 
and A describes the balance between the rate of mitosis and apoptosis (nat- 
urally cell death). The function / € C°°([0, oo)) has the following properties 

/(O) = and f'{ip) > for V > 0. (1.2) 

The tumor domain Q{t) is an unknown of the problem and, together with 
the rate at which nutrient is added to the tumor domain Q{t) and the 
pressure p inside the tumor, is to be determined. 

Three different regimes of vascularisation are introduced by the constants 
A and G : if G > and A > the tumor is low vascularised, G > and 
A <0 correspond to the moderate vascularised case, and if G < the tumor 
is highly vascularised. 

In [6] the special case / = id[o^oo) is analysed numerically. Moreover, 
in this situation, the first equation of the system is linear, and if the tumor 
domain is a sphere or an infinite cylinder, then the solution is known through 
an explicit formula. The radially symmetric case when tumors are circles is 
considered in [12], where we show that if ^ € (0,/(l)), then there exists a 
unique radially symmetric stationary solution D{0,Ra) of (1.1). The radius 
of the stationary solution depends only on A and this circular steady-state is 
exponentially stable under radially symmetric perturbations in the avascular 
case when G > 0, and unstable in the high vascularisation regime G < 0, 
result established for / = id[o,oo) ^iso in [6]. The analysis in [12] will serve 
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US in the present paper as an ancillary tool when proving the local well- 
posedness of problem (1.1) and when studying the stability properties of 
DiO,RA). 

The model, presented in [5,14,18], has been studied extensively by differ- 
ent authors, see e.g. [4,7-9,14,15] and the references therein. In particular, it 
is shown in these papers that if certain parameters belong to an appropriate 
range, then the mathematical formulation possesses a unique radially sym- 
metric solution, result matching perfectly with [12]. Moreover, the stability 
properties of this solution under general perturbations, as well as bifurca- 
tion phenomena are studied. In contrast, for the model presented in [6, 12], 
and which we consider herein, not many analytic results are available. We 
prove that also this model is locally well-posed in time, meaning that for 
appropriate smooth initial data 0,q, there exists a unique solution of (1.1), 
cf. Theorem 2.1. Though in the radially symmetric case the steady-state 
solution D(0, Ra) is exponentially stable if G > 0, we show in Theorem 2.2, 
by considering arbitrary initial data, that this solution is unstable also in 
the low vascularised case, provided G lies above a well-defined constant G*. 
This result matches the case 7 < 7* in [8, Theorem 1.2], since G is inversely 
proportional to 7. The situation when G € (0,G^=) is still an open problem. 
If G = the problem is equivalent to the Hele-Shaw problem studied in [13] 
and the exponential stability result stated by [13, Theorem 4.3] holds true. 
As a new property we establish in Theorem 2.3 exponential convergence of 
D{0, Ra) for every G > and initial data in a certain class which depends 
on G. 

The outline of the paper is as follows: we introduce in the second section a 
parametrisation for the unknown tumor domain which permites us to present 
the main results Theorems 2.1-2.3. Section 3 is dedicated to the proof of 
Theorem 2.1, and the stability results stated in Theorems 2.2 and 2.3 are 
proved in Section 4. 

2. The main results 

Let i? > be fixed for the remainder of this section. Our goal is to show 
that if the tumor is initially close to D(0,R), then problem (1.1) possesses a 
unique classical Holder solution. To this scope let h^(§), r > 0, denote the 
closure of the smooth functions G°^(S) in the Holder space G'"(S). Hereby, 
§ stands for the unit circle and we identify functions on § with 27r-periodic 
functions on JR. The small Holder spaces h^'{§) have the nice property that 
the embedding /i*" (S) is densely and compactly in h'^ (§) for all < s < r. We 
fix a G (0, 1) and we shall use functions p G V, whereby 

V := {p e /i^+"(S) : ||p||c(§) < 1/4}, 

to parametrise the boundary of the tumor domain. Obviously, V is an open 
neighbourhood of the zero function in /i^"'""(S). Given p € V, we define the 
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Figure 1. Parametrisation of the tumor domain 

C^"'"°-perturbation of the circle centred in with radius R 

Tp := {x G : \x\=R{l+p {x/\x\))] = {R{1 + p{x)) x : x G §} . 

The simply connected component of which is bounded by the curve Tp is 
the set 

np:={xeR'^ : \x\ < R{1 + p{x/\x\))} U{0}, 

with boundary dQp = Tp. Given x G Tp, the real number p(x/|x|) is the 
ratio of the signed distance from x to the circle R ■ S and R (see Figure 1). 
It is suitable to represent Tp as the 0— level set of an appropriate function. 
For this, let Np : A(3i?/4, 5R/A) ^ M be the function defined by 

Np{x) = \x\-R- Rp{x/\x\), X G A{3R/4:, 5i?/4), 

where A{3R/4:, 5i?/4) is the annulus centred in with radii 3i?/4 and 5i?/4 
A{3R/4,5R/4:) := {x G : 3R/4 < \x\ < 5R/4}. Obviously ^(3i?/4, 5i?/4) 
is an open neighbourhood of Tp and Tp = iV~^(0). Let Vp denote the out- 
ward normal at Pp. Since Tp is the 0— level set of Np, the gradient VA'^p and 
Vp must be collinear vectors. Moreover, Np is positive on the complement of 
Qp, hence Np{x + Afp(x)) > for all x G Pp and A > 0. Differentiating this 
relation with respect to A at A = yields that VNp • i/p > 0, hence 



To incorporate time let T > 0. Presuppose that the function p G C([0, T], V)n 
C"^([0, T], h^~^°'{S)) describes the evolution of the tumor, which at time t = 
is located at r2(0) = ^^p(o)- The normal velocity V{t) of the moving boundary 
Pp(t) is then given by the expression 
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This relation follows from the standard assumption that the interface moves 
along with the tumor and from relation Fp = N~^{0). With this notation 
(1.1) is equivalent to the following system of equations 



Alp = 




in rjp(t) 


Ap = 





in ^p{t) 


= 


1 


on Tp(_t) 


p = 


'-pit) 4 


on rp(j) 


dtNp = 




on rp(f) 


p{0) = 


Po 


on S, 



for all t € [0, T]. A triple {p, tp,p) is called a classical Holder solution of (2.1) 
if 

pGC([o,r],v)nci([o,r],/ii+"(s)), 

i^{-,t)p{-,t) G 6'uc2+°(17p(t)) for all t G [0,r], 

and {p,ip,p) solves the system (2.1) pointwise. Given p G V, buc^^'^{U) 
stands for the closure of BUC°°{U) in BUC^~^'^{U). Of major interest is to 
determine the mapping p which describes the evolution of the tumor. The 
functions ip and p can be then determined as solutions of Dirichlet problems, 
cf. Lemmas 3.1, 3.3, and 3.4. This is the reason why we shall also refer only 
to p as solution to (2.1). The first main result of this paper is the following 
theorem: 

Theorem 2.1 (Existence and uniqueness). Let R > 0. There exists an open 
neighbourhood O of in V such that, for any initial data G O, there 
exists a maximal existence time T := T(po) > and a unique classical 
solution p = p{-; pq) to problem (2.1) defined on [0,T(po)) which satisfies 
p([0,T(po))) C O. The mapping 

{(t,Po) : Po^O andO<t< T{po)} ^ p{t; po) G /i^+"(S) 

is smooth. 

When R = Ra, we re-discover p = 0, situation when the tumor is located 
at D{0,Ra) as the unique radially symmetric stationary solution of (2.1). 
Concerning the stability properties of this solution we already know from 
the radially symmetric case [12, Theorem 1.2] that this solution is unstable 
for G < 0. Moreover, we have: 

Theorem 2.2. Let R = Ra and > &e the constant defined by (4.21). 
Then the radially symmetric equilibrium p = is unstable for all G > G^,. 

Additionally to Theorem 2.2 we have: 
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Theorem 2.3. Let R = Ra and assume that 

1^ + ^-/(1) >0, (2.2) 
2 uo(l) 

whereby uq is the solution of (4.13) for n = 0. Given G > 0, there exists a 
positive integer Ig € N such that for alluj E (0,^o) andl > Iq we find positive 
constants Ki > and 5i > Q with the property that if \\pq\\c*+<^{S) ^ and 
Pq is 2tt/ I— periodic, then the solution p to (2.1) exists in the large, and 

IIpWIIc4+"(§) + \\p'{t)\\c^+'-{S) < ^«e"'^i/5olb4+"{s), t > 0. 
Moreover, the solution p is 2tt/ I— periodic for all t > 0. 

We will show in the Appendix that the condition (2.2) is satisfied partic- 
ularly when / = id[o^oo) ^iid = 1- 

3. The well-posedness result 

This section is dedicated to the proof of Theorem 2.1 and preparing The- 
orems 2.2 and 2.3. A fundamental difficulty in treating problem (2.1) is the 
fact that one has to work with unknown, variable domains fip. We overcome 
this difficulty by transforming problem (2.1) on the unitary disc Q := D{0, 1). 
Therefore, we define for all p € V the mapping 0p : ^ by 

Rx , / X 



Qp{x) = Rx+ -rrfil^l - 1)P rT 

where the cut-off function (/? G C°°(]R, [0, 1]) satisfies 

r 1, \r\ < 1/4, 
ip{r) = < 

[ 0, |r| > 3/4, 

and additionally max |<^'(r)| < 4. Notice that |a; — 1| > 3/4 we have Qp{x) = 
X. Given x £ S, the mapping [0, cx)) 9 r i— t- r + (p{r — l)p{x/\x\) G [0, oo) is 
strictly increasing and therefore bijective. The composition p{-/\ ■ |) has the 
same regularity properties as p on any subset of which is bounded away 
from 0, and using the chain rule we have, cf. [13], that 

K''(r))-'(r)(-^w) 

for all x 7^ 0. Consequently, Qp is a diffeomorphism mapping J7 onto fip, 
i.e. Bp e Dtff'^+'^{n, np) n Diff'^+'^{R'^,R'^). Such a diffeomorphism was first 
introduced by Hanzawa in [19] to study the Stefan problem, and it is therefore 
called Hanzawa diffeomorphism. Additionally, we have that Op(S) = Tp (see 
Figure 2). The push- forward operator induced by Qp is defined by 

@*p: BUG{Qp) ^ BUC{n), u^uoOp. 

These operators allow us to transform the problem into an abstract Cauchy 
problem over S. General results of the theory of maximal regularity, due to 
Sinestrari [23], can be used to prove existence of a unique classical solution, 
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Figure 2. The Hanzawa diffeomorphism 

corresponding to small initial data. The solution to (2.1) is then obtained 
(see Lemma 3.1 below) using the pull-back operators defined by 

where p := = (V'p, V'p)- The transformed operators A{p) and are de- 
fined as follows. Given /3 G V, A{p) : buc^~^°'{^) — )• buc°'{Q.) is the differential 
operator given by 

A{p) ■.= Q;o A oQP. (3.2) 
The operator A{p) is linear and uniformly elliptic, with 

A{p)v = bij{p)vij + bi{p)vi, Mv G 6tic^+"(0), 

whereby 

&*.(p) = V';i(ep(x))^^^i(ep(x)) + v;2(0p(^))<2(0p(^)), 

^^(/>)=V'^,ll(®p(^))+^P,22(0p(^)) 

for 1 < i,j < 2. Using (3.1) and the chain rule, we can determine the 
coefficients bij{p) and bi{p), 1 < i,J < 2, explicitly in terms of p and if, the 
cut-off function used when defining Qp. Moreover, A depends analytically 
on p 

AeC'^{V,C{buc^^''{n),buc''{n))). (3.3) 

The trace operator B : V x buc^+'^{n) x buc^+°'{n) /ii+"(S) is defined 
by the following relation 

i3{p,v,q) = Itr ^GV(e>)(e,) - v(ee(z)(ep) - ^e^, viv,(ep) 

(3.4) 

with tr the trace operator on §, i.e. tiv = v\s for v G BUC{n), and the 
curvature K^p can be expressed in terms of p by the relation 

(1 + pr + 2p'^ - (1 + p)p" 

-r.Op) = ^((l+,)2 +^,2)3/2 =• -iP)- 
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It is not difficult to see that if [p^ip^p) is a solution of (2.1) then (p,v,q) = 
(p, Qpil^, QpP) solves pointwise the following transformed problem 

A{p)v = f{v) in n, 

V = 1 on §, 

A{p)q = in n, 

( , AGR\^^ .2 ^ (3-5) 
q = k[p) — (1 + p) on S, 

dtp = B{p,v,q) on §, 

P(0) = Po. 

The notion of solution for this problem is defined analogously to that of 
solution to (2.1). In fact the problems (2.1) and (3.5) are equivalent in the 
following sense: 

Lemma 3.1. Given pQ (^V we have: 

(a) If {p^ip,p) is a classical Holder solution for (2.1), then (p, 0*^/;, 0*p) is 
a classical Holder solution for (3.5). 

(b) If {p,v,q) is a classical Holder solution for (3.5), then (p, 0^u, 6*g) is a 
classical Holder solution for (2.1). 

Proof. The proof is similar to the one in [13, Lemma 2.1]. □ 
Lemma 3.2. The mapping 

is analytic. Moreover, dK,{0)[p] = — (p + p")/R, for all p E /i^~'""(S). 

Proof. The analyticity is obvious. In order to compute the derivative one has 
only to calculate the gradient of a real valued function of three variables. □ 

We introduce now solution operators to some semilinear, respectively lin- 
ear Dirichlet problems related to our transformed problem (3.5). From the 
Leray-Schauder fixed point theorem (cf. [17, Theorem 11.3]) we obtain for 
each p G V a solution u G BUC'^~^'^{Q.p) of problem 

r Au = f{u) in Vtp, 

\ u = \ on Fp, (3-6) 

with p G V. Using the maximum principle as we did in the proof of [12, 
Theorem 2.6] we may prove the uniqueness of this solution. Consequently, 
we have: 

Lemma 3.3. Given p G V, there exists a unique solution Tip) G 6uc2+"(0) 
of the semilinear Dirichlet problem 

{A{p)v = f{v) in 0, 
1 ^ (3.7) 

V = \ on h. 

The mapping [V 9 p i— >■ T(p) G bu(?~^'^ is smooth. 
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Proof. For details we refer to the proof of [22, Theorem 4.3.5]. □ 

We consider now the solution operator corresponding to the second, linear 
Dirichlet problem in system (3.5). We state: 

Lemma 3.4. Given p £ V, there exists a unique solution S{p) € 6nc^^"(ri) 
of the Dirichlet problem 

A{p)q = in Q, 

q = K{p) ^ — on §. 

The mapping \V 3 p ^ S{p) € 5iic^"'"°(r2)] is real analytic. 
Proof. Given p € V, the mapping 

{A{p),tT)-- BUC'^+'^iVl) ^ BUC'in) X C'^+'^iE) 

is a topological isomorphism from BUC^+''{n) onto 5[/C"(f]) x C2+°(S). 
It is well-known that the function mapping a bijective bounded linear opera- 
tor onto its inverse is analytical; it can be expressed by a Neumann expansion 
in the neighbourhood of some other linear isomorphism. 

Hence, in view of Lemma 3.2 and equation (3.3) it follows that 

/ AC P^ 

S{p) = {A{p),ti)-' 0, k(p) - ^—(1 + pf 



is analytic. Since S maps smooth functions on S into BUC°°{Q) we also 
have S{p) G buc'^+°'{n) for all p e V. □ 

3.L The nonlinear Cauchy problem. We use now the solution operators de- 
fined in Lemmas 3.3 and 3.4 to transform the system (3.5) into an abstract 
Cauchy problem on the unit circle S. We put in the third equation of (3.5) 
T(p), the solution to (3.7), for v, respectively S{p), the solution to (3.8), for 
q, to obtain the following abstract Cauchy problem 

dtp = ^p), p(0) = po, (3.9) 

where 

^■):=B{-,Ti-),S{-)) (3.10) 
is a nonlinear and nonlocal operator of third order which depends smoothly 
on p. In order to prove Theorem 2.1 is suffices to show that d^{0) generates a 
strongly continuous analytic semigroup in C{h^^°'{S)) with definition domain 
/i4+"(S), that is 

-9$(0) G ?^(/i^+"(S),/ii+"(S)). 
The operator d^{0) can be decomposed as the sum of its principal part, 
which has order three in p, with an operator of first order. More exactly: 

Theorem 3.5. The operator <I> is smooth, i.e. $ G C°°(V, /i^^"(S)). Its de- 
rivative, d^{0), writes as the sum d^{0) =^1-1-^2, where 

Aip := -^d,{{A,tr)-\0,p")) for p G /i^+"(§), (3.11) 



Bi{p)v = tr (V(e»,ViVp) (Gp) and B2{p) = ^ tr (Gp, ViVp(ep)) . 
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and A2 e £(/i2+"(S),/ii+°(S)). 

In order to prove this theorem, we have to study first the regularity prop- 
erties of the operator B, defined by (3.4). It is convenient to we write this 
operator as a sum 

B{p,v,q) = ^Bi{p)v - ^Bi{p)q - B2{p), 

where Bi G {V , C{buc^+'' (Q) , h^^"^ {§)) and B2 € C"^(V, /i^+"(S)) are the 
operators defined by 

AG 
2R 

Since by the chain rule 

dmip] = - 1^1 (0)95(0) [p] - l90i(O)[p]5(O) 

+ pB^mp]no) + ^Bi(0)9r(0)[p] - dB2mp] 

for p G h^~^"{S), we must not only show that Bi, i = 1,2, have the regularity 
mentioned above but also determine their derivatives in 0. We shall see 
that the first term of this sum is the important one (corresponding to the 
operator Ai in Theorem 3.5) since it is a third order operator, and the last 
two terms are of lower order and play, as we shall see, no role when studying 
the well-posedness of the abstract evolution equation (3.9). 
Using relation (3.1) we get that 

1 + p{x) ' 

for all p € V. Particularly, we obtain find the following expression for B2 

AG AG 
B2{p) = ^ (^(1 + P)x, x) = ^{l + p), (3.13) 

wherefrom we can easily see that B2 is analytic and that 

AG 

dB2{0)[p] = ^p, ypeh'+'^iS). (3.14) 

Consider now the operator Bi. From the weak maximum principle we find 
that the function T(0) is radially symmetric, and one can easily see that 
5(0) is constant. Hence it suffices to determine dBi{0)[p]vo for a radially 
symmetric function vq E buc'^~^'^(Q) and p E h^~^°'{S). 
We state: 

Lemma 3.6. The nonlinear operator Bi is analytic, i.e. 

Bi E C"^(V X 6nc2+°(17),/i^+°(S)). 
Given vq E buc^'^" (Q) a radially symmetric function, we have that 

dBi{0)[-]vo = 0. (3.15) 



VNpiQpix)) = X - ^ ^\;^^ (-X2,a:i), x E § (3.12) 
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Proof. Let vq G buc^'^" {Q) be a radially symmetric function and p G V. In 
view of (3.12) we have that 



Slip) 



1+P 

wherefrom we obtain the regularity assumption stated in the lemma. Con- 
cering (3.15) a detailed proof can be found in [22, Lemma 4.4.2]. □ 

To finish the preparations for the proof of Theorem 3.5 one more step must 
be done. We have to determine the Frechet derivative in of the analytic 
solution operator defined in Lemma 3.4. 

Lemma 3.7. Given p £ V, the map dS{0)[p] € buc^~^°'{Q,) is the unique 
solution of the linear Dirichlet problem 

Az = in J7, 

1 AGR^\ 1 „ ^ (3.16) 

— H p p" on S. 

R 2 J ^ 

Proof. The proof is standard and we omit it. □ 

We come now to the proof of the main result of this subsection: 

Proof of Theorem 3. 5. The regularity assumption follows directly from Lem- 
mas 3.3, 3.4, 3.6, and relation (3.13). Moreover since 0i(O) = R~^dy, we 
have that 

dm[p]~,d.{{^M'\^.p")) 

(1 AG\ C AC 

+ ((A,tr)-i(o,p)) + — 5.(ar(o)M) - —p. 

The operator T defined in Lemma 3.3 can be extended to T : {p G /i2+"(S) : 
I1pIIc(§) < 1/4} ~^ buc^~^°' (Q) , because we only need there that p is of class 
C^"*"" to guarantee existence of a solution to (3.6). Whence, the operator 
defined by 

A2P := (^-^ + ^) 5. ((A,tr)"i(0,p)) + -^5.(5r(0)[p]) - ^p 

for p G /i2+"(S), belongs to £(/i2+"(S), /ii+"(S)). This completes the proof. 

□ 

We conclude this section with the proof of our first main result Theorem 
2.1. 

Proof of the Theorem 2. 1 . The key point is showing that the operator Ai , 
which can be seen as the principal part of the Frechet derivative 5$(0) 
generates a strongly continuous and analytic semigroup in C{h^~^°^{^)) for 
ah < a < 1, i.e. -Ai G ?^(/i^+"(S), /ii+°(§)). To this scope we show 
that —Ai is a Fourier multiplier. Given p G /i^"'~"(S), we consider its the 
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Fourier expansion of p = '^k&zPi^)-^'^' where p{k) := f^p{x)x^dx is the 
k-t\i Fourier coefficient of p. The well-known Poisson integral formula yields 
then 

(A, tr)-i (0, p") = -k^r\''\p{k)x^ 

for all r < 1 and a; € S. Taking the derivative with respect to r, in r = 1, we 
finally obtain 



x'' 



Y^hm^" (3.17) 

fcGZ 



for all /'(^)-^'^ ^ /i^^"(S), where or simplicity 

Afc := for keZ. (3.18) 

Following the same steps as in the proof of [13, Theorem 3.5] we obtain that 
-Ai € 7^(/i4+"(S),/ii+"(S)) for ah a e (0, 1). 

Since the constant a S (0, 1), fixed at the beginning of Section 2 was 
arbitrary we get that the assertions of Theorem 3.5 hold true with a re- 
placed by /3, for some fixed /3 € (0,a). Particularly, we have that —Ai € 
'H{h^^^{E>),h^^^{S)), for some /3 € (0, a). The definition domain of ^2 is an 
interpolation space between h^~^^{§>) and h'^~^^{S), since 

We infer from [21, Proposition 2.4.1] that the sum 5^(0) = A1+A2 also gen- 
erates a strongly continuous analytic semigroup in C{h^~^°' {§)) . In view of [1, 
Theorem 1.3.1], the set ?^(/i^+/^(S), /ii+^(S)) is open in C{h^+'^{§),h^+^{E>)), 
whence there exists an open neig hbourhood Op of in /i''+^(S) with the 
property that -d<^{p) S n{h'^+^{S),h^+^{S)) for ah p G O^. Then O := 
Ojs n h'^~^"{S) is an open neighbourhood of in V. In view of 

(/ii+^(S),/i^+^(S))e = /ii+"(S). 

we have establish that the assumptions of Theorem 8.4.1 in [21] hold and the 
proof of Theorem 2.1 follows now from this theorem. Consequently, given 
Po (z O, there exists a positive time T > and a unique classical solution p to 
problem (2.1) on [0, T] satisfying p([0,r]) C O. Moreover, the solution may 
be extended on a maximal interval [0, T(pq)) and if p is uniformly continuous 
with values in /i^+"(S), then either 

lim p(t) G do or T(po) = +00. 
t/"r{po) 

This completes the proof of Theorem 2.1. □ 
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4. Stability properties 

We study in this section the stabihty properties of the unique radiahy 
symmetric solution D{0,Ra) determined in [12] when A £ (0, /(I)) and 
G ^ 0. Therefore, we choose the constant R fixed at the beginning of Section 
2 to be i? = Rji. Particularly, functions in V parametrise domains near the 
stationary tumor -D(0, Ra)- We rediscover p = 0, situation when the tumor 
domain is the discus D{0,Ra), as the unique radially symmetric stationary 
solution of (2.1). 

In order to study the stability properties of this equilibrium we have to 
determine the spectrum of the complexification of the Frechet derivative 
9$(0), which we denote again by d^{0). The stability results established in 
Theorem 2.2 and Theorem 2.3 are then obtained by applying the principle of 
linearised stability to problem (3.9). Repeating the arguments presented in 
the proofs of Theorem 2.1 we see that the complexification of d^{0) generates 
a strongly continuous and analytic semigroup. Taking into consideration 
that the embedding /i^+°(§,C) ^ /i^+"(S,C) is compact, we deduce that 
the complexification of 0$(0) has a compact resolvent. From [20, Theorem 
III. 8. 29], we conclude that its spectrum consists only of eigenvalues of finite 
multiplicity, 

Given p G /i^'^"(S), we look in the following for the Fourier expansion of 
9<I>(0)[p]. Having shown that (9$(0) is a Fourier multiplier, the point spec- 
trum of 9<I>(0) is given by the symbol of this multiplier. The cornerstone 
of the analysis leading to Theorems 2.2 and 2.3 is the Theorem 4.1, which 
states that 9$(0) is a Fourier multiplier operator with symbol explicitly 
determined. 

Theorem 4.1. Given p G /i^"^"(S), we let p = '^k^zp{k)x^ denote its associ- 
ated Fourier series. We have that 



9$(0) 



^p{k)x'' 



.kez 



^Pkm^", (4.1) 



where the symbol {pk)kez is given by the relation 

and is the solution of (4.13) for n = \k\. Moreover, 

a(9$(0)) = {pk : keZ}. (4.3) 

In order to obtain a Fourier expansion for the derivative 5<I>(0)[p], p S V, 
we are left, cf. Theorem 3.5 and (3.17), to determine the Frechet deriva- 
tive in of the solution operator T defined in Lemma 3.3. However, the 
computations are more involved when computing dT{0), since T(0) is not a 
constant function. We have that: 
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Lemma 4.2. Given p e the function dT{0)[p] is the unique solution 

of the linear Dirichlet problem 

Aw - R\f'{vo)w = 

—2x'iLpp — 2x'^\x\ip'p + 2xiX2^p' 

- ^0,11 

2xiX2fp-2xiX2\x\(f'p+{xl-xl)ipp' 

- 2^0,12 r- To 

2x'\ipp + 2x2\x\ip' p + 2x\X2'^p' 
+ Vq,22 r- 

X\Lpp — X\ jxloj'p + 2x2'4^p' — Xl |xPw"/9 — XUpp" 

- w 

X2'^P - X2\x\lp' p -2xiLpp' - X2\x\^^" P - X2^p" . „ 

-fo,2 |^j3 m 

u; = on S, (4.4) 



where vq = T(0) and i/ie cut-off function used to define the Hanzawa 
diffeomorphism. 

Proof. The proof is standard though lengthy. For detailed calculations we 

refer to [22, Lemma 5.1.1]. □ 



The result of this lemma is not very useful yet. This is due to the fact 
that the first equation of (4.4) contains besides p and w also derivatives of 

and vq. Therefore, it is difficult to determine the Fourier expansion of 
the solution of (4.4), when knowing that of p. That is why we formally 
linearise the free boundary problem describing the stationary states of the 
full system (1.1) at the unique radially symmetric solution {Q^ipA^PA), found 
in [12, Theorem 1.1], where we simply write V'A := Tifi) oQq for the solution 
of (3.6) and pA '■= R~a ~ AGR\/4 is the composition = 5(0) o 0q. 
By doing this we shall find a nice decomposition of the derivative 9T(0) as 
a sum of two operators (see Lemma 4.3 below). Hence, we consider now 
perturbations of the radially symmetric solution of the form 

Ipe = 1pA + £1p, 
Pe = PA + ep, 

= {re'' : < r < RAil + ep{s)), s gR} , 

where we simply write p{s) = p{e'') for all s G M. Here, e is a small param- 
eter, and tp,p,p are unknown functions. The linearisation of problem (1.1) 



WELL-POSEDNESS AND STABILITY PROPERTIES 



15 



is then the fohowing free boundary problem 

AV' = /'(Va)-V' mD{0,RA), 

Ap = mD{0,RA), 

^ = -RaiP'a{Ra)p ondD{0,RA), 

P = -(^ + ^^)p-^P" ondD{0,RA), 



Ra 2 J' Ra 

Gdytj) - dyp = ^^^^ p on dD{0,RA)- 

(4.5) 

We look now for a connection between the linearisation (4.5) and the Frechet 
derivative of $ in 0. To this scope, we transform first the system (4.5) to 
the unitary disc, i.e. we set 

w{x) = iP{Rax), z{x) = p{Rax), vo{x) = TpAiRAx), 

for X € 0, and substitute these expressions in (4.5). This leads us to the 
following system of equations 

Aw = R\f'{vQ) -w in n, 

w = —dyVQp on S, 

Az = Q in 0, 

(4.6) 

/ 1 AGR\\ 1 „ 
G 1 AGRa „ 

^ ^^^^ - lu^-'P = °" ^- 

Given p G /i^+"(§), we let W{p) G buc^+°'{n) denote the solution to the 
linear Dirichlet problem 

( Aw = R\nvo) -w in fi, 
1 w; = —d^vop on S. 

Further on, we want to determine a relation between dT{0)[p] and W(p). 
Therefore, we define the extension operator £ : h^~^"{S) — t- 6nc^"^"(r2) by 

£{p){x) := L i(^)f^ + vo,2{x)^) P (^) V{\x\ - 1) (4.8) 

V 1^1 fI/ VfI/ 

for X € ri. Using these operators we can now write dT{0) as the sum of 
W and £. This decomposition is very useful because we got rid, in this 
way, of all the terms from the right hand side of first equation of (4.4). 
Recall that our goal is to determine the Fourier expansion of du{dT{0)[p]) 
when p G h^~^'^(S). It turns out, cf. Lemma 4.3, that du{£{p)) is collinear 
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with p for all p G /i^''"°(S). Moreover, using ODE-tecliniques we are able to 
determine an expansion for W(/9) for all p G cf. (4.15). Indeed, we 

have: 

Lemma 4.3. It holds that dT{0) = W + Moreover, given p G /i^"'"°(§), we 
have that 

d,{£{p)) = aAP, (4.9) 

where oa '■= d'^U/dr'^{l,R'j^) > and U is the solution of the parameter- 
dependent problem: 



^ ^(r X) + - — (r X) 
C/(1,A) 



A/([/(r,A)), 0<r<l, 

0, 
1, 



(4.10) 



with A G [0, oo). 



Proof. The proof follows by direct computation by taking into consideration 
that Wo = r(0) = C/(| • |,i2^). □ 

In virtue of Lemma 4.3, if we determine a Fourier expansion for W(/9), 
then we completed the task of determining the expansion of 9T(0) [p] for all 
p G h'^'^^i^). For this reasoning, we consider expansions of the form 

W(rx) - 4. ^^^^fc 
p{x) 



l.lll 



with r G [0,1] and x G S. Substituting these expressions into (4.7), and 
comparing the coefficients of , we come to the following problem for the 
unknown function Aj. : 



1 



r 

Akil 



(4.12) 



R\f'{vQ)Ak, < r < 1, 
-v'o[l)m. 

We have used here the relation d^VQ = v'q{1) on S, where we identified vq 
with its restriction to the interval [0, 1]. In order to prove the existence and 
uniqueness of the solution to (4.12) we consider first the associated problem 
(4.13). The solution of (4.12) will be then expressed in terms of the solution 
of this new system. Thus, given n G N, the problem 



,, 2n + l , 
u" H u' 



u(0) 
n'(0) 

has a unique solution Un G C°°( 



= 1, 

= 0, 



< r < 1, 



1.13) 
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With this notation we have: 

Lemma 4.4. Given A; € Z, problem (4.12) possesses a unique solution Ak G 
C°°([0, 1]) explicitly given by 



Ak{r) = ^^p{ky^mi^)^ 0<r<l. (4.14) 
'^\k\{i') 

We denoted by Un,n € N, the solution of (4.13). 



Proof. That is a solution of (4.12) fohows by direct computation. The 
uniqueness may be obtain by a contradiction argument. □ 

We give now a short proof of the theorem stated at the beginning of this 
section. 

Proof of Theorem ^.1. From Lemma 4.4 we obtain the following expansion 
for W{p) 



W{p){rx^) = Y, 



u\k\{r) 



p{k)x^ (4.15) 



for all p € /i^"'"°(S). Let us determine now the constant a a form Lemma 
4.3. From the first equation of (4.10) we find that the constant satisfies 
the relation aA = I^Afi^) ~ ^o(-'-)- Moreover, from the same equation and 
relation (3.27) in [12] we have that 

v'^{l) = AR\/2, (4.16) 

hence 

aA = R\ (/(I) - ^) . (4.17) 

In view of Theorem 3.5, 3.17, and (4.15) we conclude (4.1) and the proof is 
complete. □ 

4.1. Estimates for the symbol of the derivative of In order to study the 
stability properties of the unique radially symmetric equilibrium determined 
in [12] we need to study the sign of symbol {pk)k&z given by (4.2) in depen- 
dence of the parameters (A, G) € (0, /(I)) x (0, oo). We consider in here just 
the case when G > 0, since for G < we already established in [12, Theo- 
rem 1.2. (d)] that this circular equilibrium is unstable. It is worth noticing 
that //fc = p_k for all A; € Z, so that we consider just the terms p^ with 
A; G N. At first we have to as certain that is in the spectrum of (9<I>(0) for 
all (AG) G (0,/(l)) X [0,oo). 

Proposition 4.5. We have that 

m = 0. 

Proof. Since vq = T(0) satisfies vo{r) = U{r,E?)^ for all < r < 1, we have 
that 

v'^ + \'^ = R\f{v^) in (0,1). 
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Differentiating this equation with respect to r and setting wq '■= v'q we get 



1 



R\f'{vo)wo, 



(4.18) 



[ wo{l) = v'oil). 

Given r G [0,1], define w{r) = VQ(l)rui{r)/ui{l), where ui denotes the 
solution of the problem (4.13) when n = 1. With c := 'yQ(l)/ni(l), we obtain 
by differentiation that 

w'{r) = c{ui{r) + rui{r)), w"{r) = c{ru'{{r) + 2n'^(r)), 

which in turn imphes that w is a solution of the equation 



.^1 / 

w H — w 



1 



-w = c{ru'[{r) + 3n'i(r)) = Rj^f'{vo)w. 

Moreover, for r = 1, we get w{l) = "^0(1). Since the solution of (4.18) is 
unique, we get wq = w, meaning that Wo(r) = VQ{l)rui{r)/ui{l) for all 
< r < 1. Rearranging this relation it follows that 

v'oir) _ ui{r) 

which leads to 

u[{r) _ v'^ir)r v'^il) - v'^ir) v'q{1) _ 



1 



ni(l) v',{lfr^ v',{l) 

If we choose r = 1 this last relation leads to 

ni(i) - .^(1) - .^(1) 

In view of relations (4.16) and (4.17), and the definition of oa, we have 



1. 



AR\ 



-2 + 



2/(1) 
A 



Inserting this result in the expression of ^ui yields 

^A 



/^i 



1 1 

r^^r:^ 



G 



A 



2 + ^)+^-/(l) 



0. 



□ 



Proposition 4.5 reveals that /ii = belongs to the spectrum of 9<I>(0). 
However, we show now that the sequence fik —^k^oo —00. This is obviously 
true if the sequence (n^l)/ii„(l))„gfs} is bounded. Even more, we have: 



Lemma 4.6. R holds that 



Unil) 
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Proof. Let n G N and set v := Un+i — Un ^ C°°([0, 1]). Recall that Un is an 
increasing function for all n € N. This can be seen since 

n^{r) = l + 1^ ^ l\'^+^nvo{T))un{T)dTds, 0<r<l, (4.19) 

relation which is obtained by multiplying the first equation of problem (4.13) 
with r^""''^ and then integrating twice. From (4.13) we obtain 

9r) -I- 1 9 
v" + =^^^v' = Ry'ivo)v--u'^^^, 0<r<l, 

V(0) = 0, 

v{0) = 0. 

Furthermore, we have 

lim^ = lim<+^^^^-<^^^=i^i/(.o(0))^ 1 1 



t^o t t^o t V2n + 4 2n + 2 

= -^^^''''°'°» (2„ + 4H2n + 2) < 

which implies by (1.2) that v'{t) < for t G (0, 6) and some 5 < 1. Thus, v is 
decreasing on (0, 5). Let now t € [0, 1] and mt := minjg f < 0. A maximum 
principle argument shows that the non-positive minimum must be achieved 
at t, mt = v{t) which implies Un+i{t) < Un{t) for all t G [0, 1]. Particularly, 
Un+i{l) < iin(l)- It also holds that 



I " /T\ 2n+l 



dr 



>Rl (-) f'{vo{T))Un+l{T)dT 

for all t G [0, 1], hence u^_|_;^(l) < ti(^(l). Moreover, from 

/*T2"+V'(t;o(r))n„,(r)dr 
JO 



<(t) 
it follows that 



^2n+l 



<(l)=iii j\^-^^f'{v,{T))Un{r)dT. 

Jo 

The dominated convergence theorem implies n^(l) \ 0. A similar argument 
provides ^^(1) \ 1. Therefore, n^(l)/ti„(l) ^0, and we are done. □ 
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4.2. Instability for G > G*. We prove now that if G is large enough, then the 
radially symmetric equilibrium D(0, Ra), determined in [12, Theorem 1.1], is 
unstable, in the sense that problem (3.9) has nontrivial backwards solutions. 
This reveals that the exponential stability result stated in [12, Theorem 1.2. 
(a)] gives a false impression about the stability properties of the radially 
symmetric equilibrium. If A; € N is large enough, we find in view of the 
Lemma 4.6 a unique value of the parameter G such that /i^ = iff 
G = Gk- More precisely, given A; € N such that 

4^ + ^-/(1)7^0, 

2 Ufc(l) 

we set 

LfcS + J-k 

G, := (4.20) 

2n,(l)+^ ^^^^ 
To prove the statement of Theorem 2.2, we set 

:= min {g,. : + ^ " /(I) < o| > 0. (4.21) 

I 2 Uk{l) J 

Moreover, since Gk -^k->oo oo, the minimum must be achieved, i.e. we find 
at least an integer such that 



^ <o(l) 
2«fco(l) 



+ ^-/(l)<0, (4.22) 



and G^: = Gk^- With these preparations we are able to prove the instability 
result. 

Proof of Theorem 2. 2. Let G > G* be given and be an integer such that 
G* = Gfcg and (4.22) holds true. It follows that 



which implies > 0. We are left to check the following instability assump- 
tions 

J CJ+ (a$(0)) =a(9$(0))n {AG C:ReA>0} / 0, 

\ inf{ReA : A G (j+(5<5(0))} > 0, 

where d^{0) stands here for the complexification of d^{0). The first one 
is clear since /ifc^ > 0. Moreover, we infer from Lemma 4.6 that /Xfc — 7'fc_j.oo 
— oo, and therefore the unstable spectrum (T+(5$(0)) contains finitely many 
positive eigenvalues of the Frechet derivative d^{0). Thus, we found out 
that the assumptions of [21, Theorem 9.1.3] are satisfied and therewith the 
radially symmetric solution D{0,Ra) is unstable. □ 
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4.3. Exponential convergence of 27r/l— periodic data. In this section we show 
that the exponential stability result stated in [22, Theorem 3.0.3 (a)] can be 
generalised for solutions of the original problem (1.1) which correspond to 
initial data pQ closed to the unique radially symmetric equilibrium and 2tt/1- 
periodic. The positive integer I depends on the constant G, which is chosen 
to be positive, so that D(0, Ra) is the unique equilibrium of the problem 
(1.1). The main result of this section, Theorem 2.3 requires the following 
assumption 

which means that the eigenvalue fiQ = /io(G) is negative for all G > 0. This 
assumption is satisfied for example if Ra = 1 and / = id[o,oo) (see Appendix). 

In order to prove Theorem 2.3 we introduce first appropriate subspaces 
of the small Holder spaces. Given € N and Z S N, n > 2, we define the 
subspace of /i'^'"^°(S) consisting of 27r/Z-periodic functions by 

/if+"(S) := {p G /i^+"(S) : p{x) = pie^-^'l^x) for all x € S}. 

Set further := V n /if+°(S). The Fourier series associated to p G /if+"(S) 



IS 



p = Y,p(ki)^''\ 

k=0 

where p{kl) is the kl—th. Fourier coefficient of p. Our first objective is to 
prove that if /> G V/, then $(p) G /i;^+"(S), that is 

G C7~(Vi,/i^^+"(S)). (4.23) 

Having shown (4.23), by choosing / large enough we can exclude the eigenval- 
ues /ifc with k small from the spectrum of (9<I>(0). These are the eigenvalues 
which we could not estimate whether they are negative or not. In this way 
we also eliminate pi = from the spectrum. 

Let /? G V/ be given and let ip := Q^v, where v is the solution of (3.7). We 
prove that ip satisfies tl'{x) = 0(e^'^*/'x) for all x G ilp. Therefore, we must 
prove first that, if x G Op, then xe^'^*/^ belongs to Op. Indeed, given x G Op, 
we have that 



xe^^^/'l = \x\<R[l + p{ — ]]=R\l + p 



2-Ki/l 



x\) ) \ \ Ixe^'^*/' I / / ' 

which implies that xe^'^*/' G Op. Recall that the function ■0 G 6uc^'^"(0p) is 
the unique solution of the Dirichlet problem 

f AVj = /(Vj) in Op, 
I V = 1 on Fp. 



22 J. ESCHER AND A.-V. MATIOC 

Setting tlj{x) := 'i/j(xe^'^*/') for x G Op, we have defined in this way a further 
solution of (4.24), since 

A^(x) = AV^(xe2-/') = /(^(xe^-/')) = f{^{x)), 

for all X € Op and ip = 1 on Tp. The uniqueness of the solution to (4.24) 
implies that if) = ipi. Thus ^(x) = ipixe'^'^^^^) for all x € Op. Following 
the same schema we can prove that p := 0*q, where q is the solution of 
(3.8), satisfies p{x) = p(e^'^*/'x) for all x E Op provided p € V;. With these 
preparations we state: 

Lemma 4.7. Given I > 2, the operator $ maps smoothly Vi into hj~^°'{S). 

Proof. It remains to show that Bi{p)T{p) and Bi{p)S{p) are 27r//-periodic. 
The assertion follows then in view of (3.10) and relation (3.13). We prove 
just the assertion for Bi{p)T{p), the proof that Bi{p,S{p)) is 27r/Z-periodic 
follows analogously. Indeed, given x G S we have 

fii(p,r(p))(x)=(v(e»(ep(x)),viVp(ep(x))) 
=a,^(Gp(x))|viVp(ep(x))|, 

where ^ is the solution of (4.24). Additionally, | VA'p (Gp (xe^''*/')) | = 

\S/ Np(@p{x))\ for all x E S. In order to prove that also 9^^/;(Gp(xe2'^*/')) = 
d^tp^Qp^x)), X € S, we introduce the rotation matrix 

cos(^) -sin(2f) X 

V -n(f) cos(^) )' 

We infer, from ^(x) = V {xe^'^^/^) = V'(M-x), that Vi{j{x) = • V'(/'(M-x) 
for all X € Op. Thus, 

(Op (xe^-/')) = (VV(ep(M.x)),i/p(ep(M-x))) 

X - p'(x)(-X2,xi) \ 



Vij{M -Qpix)) ,M ■ 
M ■Vtl;{Qp{x)),M ■ 
V'0(ep(x)),M^ -M 



^(l + p)2+p'2(^) 
X - p'(x)(-X2,Xl) \ 



^{l+p)^+p'2(x)/ 
X - p'(x)(-X2,Xi) 



V(l + p)2+p'2(x) 

(v^(ep(x)),i.p(ep (x))) 
(Gp (x)) 
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for all X eS. Summarising, Bi{p,T{p)) {xe^'^'l^) = 6i(p,r(p)) (x) for all 
X G S, and the proof is completed. 



□ 



The Frechet derivative (9$(0) of the mapping IS, m 

view of (4.1), given by the relation 



Vk&L 



x^\ 



where {fikl)k are defined by (4.2). We come now to the proof of the expo- 
nential stability result for 2tt/1— periodic data: 

Proof of Theorem 2.3. Let G > be given. Since p^i^^ — ^|fc|->oo —oo, we 
find a positive integer Iq such that < fiQ for all |A;| > Iq- Let / > Iq 
be fixed. Li view of relation Lemma 4.7, we find that the restriction $ € 
C°°{Vi,hl+"{§)) satisfies the assumptions of [21, Theorem 9.1.2]. Indeed, 
since / > it holds that pkl ^ Po for all € N. Consequently, the spectrum 
of the complexification of the d^{0) consists only of the negative eigenvalues 
{pkl ■ k G N}, and is bounded away from the positive half plane by /ig- The 
assertion follows now immediately from [21, Theorem 9.1.1]. □ 

5. Appendix 

We show now that the condition (2.2), meaning that Pq{G) < for all 
G > 0, is not to restrictive. 



Observation 5.1. The assertion 

Au',{l) 



+ A-fil)>0 



2 no(l) 

is fulfilled when f = idjg.oo) ^'^'^ = 1- 

Proof. In view of Proposition 4.5, our assertion is equivalent with 

+A- /(I) >0 = p,{G) = 4^ + A- /(I). 
2 uo{l) 2 ui{l) 

Consequently, we have to show only that 

^o(l) ^ <(1)^ 

Mo(l) ""1(1)' 

We assume now that uq, the solution of (4.13) when n = 0, is analytic 
and the Taylor series associated to uq in 



^0 = Qfc 



k=0 



24 J. ESCHER AND A.-V. MATIOC 

converges on [0, 1]. Problem (4.13) writes now as follows 
XUq + u'q — XUq = 0, < x < 1, 
no(0) = 1, 
n'o(O) = 0. 

From the initial conditions of (4.13) it follows immediately that uq = tto(O) = 
1 and oi = Uq{0) = 0. Plugging uq and its derivatives in the first equation 
of the system, one finds out that 

OO CXD OO 

k{k + l)afc+ix'= + ^(/c + l)afc+ix'= - aox - ^ a^-ix'' = 0. (5.1) 

fc=l k=0 k=2 

Identifying the coefficient of in (5.1) yields 

ao 1 

= T = 4' 

[k + ly 

and, from oi = 0, we deduce that 

k ^ 

a2k+i = and a2k = 7^^' ^ ^' 

n=l 

thus, 

uo(x) = l + J](n (^j^'', xe[0,l]. (5.2) 



We make now the same assumption on m, the solution of (4.13) when 
n = 1. We then get, that ui is the solution of the following system 

xui + 3u'i — xui = 0, < X < 1, 

ui(0) = 1, 

ii;(o) = 0. 

As above, we obtain 

OO OO OO 

^ k{k + l)ak+ix^ + 3 ^(/c + l)ak+ix^ - ao^ - ^ ak-ix^ = 0, 

k=l k=0 k=2 

and therefore 

k ^ 

a2k+i = and agfc = TT — — — — , V/c G N. 

Iniln + 2) 

n=l 

Thus, 

OD / k ^ \ 
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It is worth noticing that the Taylor series associated to uq and ui, respec- 
tively, in define analytic functions on the whole real line, so that the rep- 
resentations (5.2) and (5.3) are valid. 
With three exact decimals we have that 

^ « 0.446 > 0.240 
which leads to the desired conclusion. □ 
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